Investigation of Possible Dark Matter Direct Detection in Electron Accelerators 
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We investigate a possibility of neutralino dark matter (DM) direct detection in the future electron 
accelerators. That is counting of high pr electron recoil events by neutralinos in halo. If selectron 
and neutralino masses would be precisely measured in future collider experiments, the beam energy 
could be tuned so that the scatterings are dominated by on-pole selectron exchange. When selectron 
and neutralino mass difference is smaller than O(10) GeV, the elastic cross section exceeds over 
micro barn. Discovery of the high pr electron events would be a firm prove of the neutralino DM 
component in halo. In the experiment, the electron beam energy must be tuned within O(10) MeV 
and the electron beam with high currents of O(100)/l is required for the detectors of the total length 
of a few hundred meters so that the sufficient event rate is obtained. The dependence of the event 
rate on the DM velocity distribution in halo is also discussed. This method might be applicable to 
other DM candidates. 

PACS numbers: 95.35.+d, 95.55.Ka, 11.30.Pb 



Nature of the dark matter (DM) in the universe is an 
important problem in particle physics, astronomy and 
cosmology. The lightest neutralino, x , m the minimal 
supersymmetric (SUSY) extension of the standard model 
(MSSM), is a good DM candidate since the lightest SUSY 
particle (LSP) is stable due to the R parity The cos- 
mological DM abundance, which is now precisely mea- 
sured by the WMAP 2] , constrains properties of the neu- 
tralino and the SUSY particle mass spectrum if the DM 
is generated in the hot thermal bath in the early universe 
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The neutralino works well as the cold DM in the struc- 
ture formation in the universe. High resolution iV-body 
simulations show that the cold DM hypothesis explains 
well the large structure of our universe [j] . On the other 
hand, the DM distributions in smaller scales than O(l) 
Mpc are still unresolved. The local DM abundance and 
the DM velocity distribution on the neighborhood of the 
solar system are also not well constrained from the rota- 
tion curve measurements. The direct DM detection on 
the earth and the indirect detection of anomalous cosmic 
rays produced by the DM annihilation may give clues to 
the problems. 

The conventional DM direct detection relies on nuclear 
recoil in nuclei-DM elastic scattering. The sensitivities of 
the proposed experiments cover a significant part of the 
MSSM parameter space. However, to evaluate the local 
DM density and velocity distribution from the counting 
rates, precise determination of the hadronic matrix ele- 
ments is necessary. The cross section also depends on 
various SUSY parameters such as relatively suppressed 
Yukawa coupling of strange quark (or tan 0) , heavy Higgs 
mass and Higgs-neutralino-neutralino coupling and so on. 
To this end, positive signals in the direct detections may 
not necessarily prove whether the DM is the neutralino 
since they merely measure the DM-nuclei scattering cross 
section. 

The nature of the LSP would be measured if the SUSY 



particles are discovered in future high energy colliders. 
This is because the LSP will be copiously produced 
from the cascade decays of the heavier SUSY particles. 
The LHC experiment is scheduled to start on 2007, and 
squarks and gluino with masses up to 2.5 TeV can be 
discovered. Furthermore, the interaction of the lightest 
neutralino would be also measured at the LHC Q and 
a future linear collider (LC) 0. While they might be 
successful to determine the thermal relic density of the 
universe and provide consistency checks of the neutralino 
DM assumption, it relies on the assumptions that the 
neutralino is stable in the cosmological time scale and 
the thermal history of the universe follows the standard 
big bang scenario. It is important to observe the neu- 
tralino DM in more direct and less ambiguous ways. 

In this paper we investigate a possibility of direct neu- 
tralino DM detection which might be realized in the fu- 
ture electron accelerators. The electrons in the beam can 
interact with the DM neutralinos in our neighborhood. 
The electron-neutralino elastic scattering is induced by 
the s-channel exchange of selectrons e - , which are su- 
perpartners for electron, in addition to the i-channel Z 
gauge boson exchange. If the beam energy would be 
tuned to the difference between the neutralino and se- 
lectron masses (Am), the elastic scatterings between the 
DM neutralinos and the beam electrons are dominated by 
on-pole selectron exchange and the cross section is sup- 
pressed only by square of Am. Building of such an elec- 
tron beam might be considered when the mass difference 
is precisely measured in the future collider experiments 
such as LC experiments. 

If a high intensity electron beam would be available, 
one could determine the nature of the DM without any 
ambiguities by detecting on-pole production of selec- 
trons, because the measurement proves that the DM 
is the neutralino LSP. Note that the relevant couplings 
would be directly measured at LC when selectron pro- 
ductions are accessible. Therefore, the DM physics that 
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might be explored by a high intensity electron beam 
is unique. However, one needs several technical break- 
throughs to realize it. We clarify the requirements to 
the electron beam and detectors in this paper. It is also 
shown that the number of events might have a sensitivity 
to the DM velocities when the selectron decay width is 
suppressed by the mass difference or the coupling. This 
implies that if there would be enough statistics, the DM 
velocity distribution in the halo and the DM wind due 
to a circular motion of the solar system might be con- 
strained by this experiment. 

First, let us discuss the elastic scattering process be- 
tween electrons in the beam and the DM neutralinos 
and evaluate the expected number of events. The DM 
neutralinos are highly non-relativistic in the universe, 
and the local DM velocity in our neighborhood is typ- 
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by the on-pole selectron exchange. In this case the spin- 
averaged differential cross section with respect to the an- 
gle between the beam and the recoiled electron, 6, is given 
as 
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with s square of the center-of-mass energy and Tg- the 
selectron decay width. Here, the function ^(cos^) is de- 
fined as 
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A(cqhB) = — x_(l-cos0), (2) 
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and the energy of the recoiled electron, i£ r ecoiii is also 
given by it as 



^recoil = #beam(l + A(cOS0)) 1 . 

When x° is bino-like, the selectron decay width is 
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where Y is the hypercharge for e~ and On is the neu- 
tralino mixing matrix element. For simplicity, we take 
Ow — \ and Y = — 1 for the right-handed selectron in 
the following. In this case 
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When the neutralino and selectron masses are close to 
each other, we find from Eq. that the cross section 
at E heam = £bcam(- Am) is suppressed only by (Am) 2 . 
If the DM velocity dependence of the cross section is 



negligible, the expected number of events, N, is 
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Here, j is the beam current, L the detector length, and 
T the duration of experiment. We fix the local DM mass 
density, pdm, to be 0.3GeV/cm 3 in this paper. If the 
beam can be polarized, the expected number of events is 
multiplied by a factor two. This elastic scattering process 
is also noticed in an earlier work in Ref. 0. 

When Am/m^.o 1, the elastic cross section can be 
enhanced if the beam energy could be tuned to -Ebeam, 
however, the requirements for the DM direct detection 
in electron accelerators would be severe. We now discuss 
the requirements to observe such events in order. 

First, for our purpose, the neutralino and selectron 
masses and the coupling have to be measured with suffi- 
cient precision at the earlier collider experiments so that 
-Ebeam and Tg- are determined. Especially, the uncer- 
tainty of -Ebcam must be smaller than r§- . Otherwise, the 
beam energy could not be tuned, so that the elastic pro- 
cess between the DM neutralinos and the beam electrons 
is not dominated by the on-pole selectron exchange. This 
implies from Eq. (J5J that -Ebeam should be determined at 
least with precision of O(10 -3 ) for Am/m^o = 10%. 

At the LHC, the mass difference between the neu- 
tralino and sleptons might be measured with the error 
on the order of a few GeV by using the events with 
lepton-antilepton pair for favorable parameters @. The 
LC would be able to measure the absolute LSP and slep- 
ton masses with the error O(50) MeV using the threshold 
scan and the end point measurements 9]. The error for 
the mass difference Am, determined by the end point 
measurement, may be even smaller. When the beam en- 
ergy spread is a dominant uncertainty in determination 
of Am, Am might be determined with the precision of 
10 -3 at the LC. The neutralino interaction would be also 
measured precisely at the LC so that the selectron decay 
width might be determined. 

Next is the requirements for the electron beam. The 
energy spread of the electron beam must be less than 
O(10 -3 ) for Am/m^o = 10% by the same reason as for 
the determination of -Ebcam- In addition, very high cur- 
rents are required in our proposal as seen in Eq. JHJ), 
which is about 10 times higher than that for the cur- 
rently proposed Super B factory. 

In the KEKB at KEK, which is an asymmetric 
electron-positron collider for B physics, the averaged 
beam currents in the low energy positron (3.5 GeV) and 
high energy electron rings (8.0 GeV) are 1.861A and 
1.275,4, respectively. The PEP-II at SLAC also has com- 
parable currents. Now upgrade of the KEKB to the Su- 
pcrKEKB is proposed, aiming for the beam currents 9.4 
and 4.1,4 [lj|. While the energy spreads of these acceler- 
ators are smaller than 10~ 3 , their beam currents are not 



3 



much enough for our purpose. 

The beam focusing at the interaction point is the 
sources of the beam instability in collider experiments 
with high luminosities. For the DM direct detection in 
electron accelerators, it would not be a serious problem 
since the beam focusing is not required. Rather, the syn- 
chrotron radiation (SR) from the electron beams would 
bound the beam currents. The SR at the arc sections 
of the accelerators may damage the beam pipe and also 
causes the beam power loss. 

A hint to solve this problem may be in a technology 
called Energy Recovery Linac (ERL) In this scheme, 
the electron beam energy is lowered by transferring the 
energy to the RF power, and the power is used to ac- 
celerating the electron beam again. The principle of the 
ERL technology have been tested at Thomas Jefferson 
National Accelerator Facility ^2] ■ Various facilities using 
the ERL technology, such as photon factories, electron 
cooling and so on, are proposed. If it would be possible 
to have a storage ring where the beam energy is lowered 
before the arc section to keep the beam power in the RF 
for accelerating the electrons in the straight section, one 
could achieve the high current electron beam with less 
power consumption [20|. 

Even if an electron beam of O(100)A is possible, the 
number of events would be only around 1 event/lOm/year 
for mass difference of the order of 10 GeV. Note that 
because this is a fixed target experiment where the target 
is on the beam line, many detector units should be placed 
along the beam line. The total length must be as long 
as 300 m— 1 km. This means that each detector units 
should be as simple as possible. 

We note that the signal events should be distinguished 
over the backgrounds (BGs). When Am/m^a <C 1, the 
recoiled electrons in the signal events have energies be- 
tween i?bcam(l — 2Am/m^o) and -Ebeam and large trans- 
verse momentums. This is because A(cos 9) in Eq. (J2J) 
is suppressed by Am/m^o and the angular distribution 
of the signal electron is almost spherical. Furthermore, 
the momentums of the signal electron must be pointing 
to the beam line. The kinematics of the signal events is 
well-constrained. 

The signal electrons would be detected by placing 
either electromagnetic calorimeter or tracking chamber 
with solenoid magnets along the beam line. The expected 
BGs would come from either cosmic rays or the beam in- 
teraction. To remove accidental cosmic ray BGs, we need 
a reasonable resolution for electron tracks to reject any 
tracks crossing the beam line from outside of the detec- 
tor. 

Other BGs are electrons scattered by the beam gas. 
They could be reduced by measuring the electron mo- 
mentums since the BG electrons are forward-going and 
Pt <C fibcam- The electron momentum perpendicular 
to the beam axis can be measured precisely with mag- 
netic fields, while the resolution of the electron momen- 
tum along the beam axis would be much worse than that 
in ordinary wire chambers. 



Since the magnetic field along a long detector will 
destabilize the beam, the detector should be divided to 
short segments along the beam line so that the direction 
of the magnetic field in the each segment is reversed to 
the next segment. This segmentation would also suppress 
the BG electrons due to pile up from the upper currents, 
and the background level might be lowered. 

7r~s produced from photo- nucleon interactions in the 
beam pipe would be also a source of the BGs. They 
are even more low energy compared to the electrons 
from beam-beam gas interactions, and multiply pro- 
duced. The detectors with tt~~€~ separation abilities 
may be useful to reduce the BGs furthermore. Transi- 
tion radiation detectors (TRD) may have a benefit for the 
purpose. It is composed of layers of radiators and X-ray 
detectors. The radiators emit soft X rays when charged 
particles are injected. The radiation efficiency for tt~ 
is lower than that for e~ for 0.5 GcV< p < 100 GeV, 
and the information may be used to discriminate 7r~s 
over e~s. The rejection efficiency of about 99% may be 
achieved for a total TRD width ~ 50cm Q . 

We have sketched a rough design for the detector. 
While the signal events have the well-defined kinemat- 
ics and they might be discriminated from the BGs by 
such a detector, the signal event rate would be small and 
we need a long duration of the experiment. The per- 
formance of the detector must be studied very carefully, 
because the dominant part of the BGs must come from 
mis- measurements. 

We have neglected the dependence of the DM veloc- 
ity distribution in our neighborhood on the event rate in 
the above discussion. However, when the selectron de- 
cay width is small, the dependence cannot be neglected. 
When -Ebcam = i^bcam, the deviation of y/s from mg- in 
the signal process is typically 

~ (i^v) m 

with (vn) the average of DM velocity along the beam 
axis. It is found from Eq. Q that this value is com- 
parable to or larger than the selectron decay width for 
Am/m^o^ 10%. Thus, the DM velocity distribution 
affects the expected event rate when Am/m^o < 10%. 

If the DM neutralinos are discovered in the electron 
accelerator experiment, we might constrain the DM ve- 
locity distribution by collecting the signal events. Since 
our DM direct detection in the electron accelerator needs 
several breakthrough technologies, it might be premature 
and speculative to discuss the measurement of the DM 
velocity distribution in addition to the DM detection. 
However, since technologies to measure the DM velocity 
distribution in the conventional DM detections are still 
limited, especially in cases of the small counting rate |l5j . 
the measurement of the DM velocity distribution using 
the electron accelerator might be considered as one of the 
alternative possibilities. 

In the following we evaluate the averaged cross section 
assuming the spherically symmetric isothermal sphere 
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(SSIS) model for the DM velocity distribution. The DM 
velocity distribution in the rest frame of the earth is then 
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We take the velocity dispersion of our local halo as 
an = 270km/s. v e h is the velocity of the earth with re- 
spect to the halo, and it is given as v e h — v es + v s h, 
with v es and v s h the velocities of the earth with re- 
spect to the sun and of the sun with respect to the 
halo, respectively. The solar system is moving toward 
the constellation Cygnus ( (a, <5)=(21M2m01.053sR.A., 
+48°19'46.71"decl.) (J2000.0) in the equatorial coordi- 
nate system) with speed \v s h\ = 233km/s. The earth's 
speed in the orbital motion is \v es | = 29.8km/s 16]. The 
motion of the solar system generates the DM wind from 
the constellation Cygnus to the observers on the earth. 
When the beam axis is parallel (perpendicular) to the 
DM wind, («„) is y/°i/3 + vi h {yffifi). ^ 

Being the experiments on the ground, the beam axis 
rotates around the the earth rotation axis. It leads to 
modulation of the event rate with a period a sidereal day 
(23/i56m4.09s) due to the DM wind. Let us assume that 
the azimuth angle from south for the beam axis is and 
the experimental site is placed at latitude ip. The angle 
between the beam axis and the direction of the DM wind, 
0, is modulated as 

= cos <5* cos 8 cos(t — a — t*) + sin <5* sin. 6 (9) 

where sin 8* = cos 9 cos tp and sin f * = — sin 9 / cos 8+ , and 
t is the local sidereal time at the experimental site. 

In Fig. ^ the modulation of the cross section for the 
signal event during a sidereal day is presented in cases of 
8* = 0° and 45°. Here, = lOOGeV, Am = lOGeV, 
and E'beam = -Eboam- The cross section averaged by 
Eq. (jSJ is reduced from that for the on-pole selectron 
exchange process (~ 12/^barn). The latitude for the ex- 
periment site tj) is assumed to be 0° for simplicity. The 
phases for 8+ = 0° and 45° are not equal to each others 
when ijj is different from 0°. (See Eq. ©.) Since 8 is 
close to 7r/4, the amplitude of modulation in 8* = 45° is 
almost maximum. The maximum (minimum) point on 
the curve corresponds to a case that the beam is almost 
perpendicular (parallel) to the DM wind. The angle a 
can be inferred by the phase of the observed modulation. 
It might be difficult to constrain 8 using a single beam 
line with a fixed energy due to the parameter degeneracy 
with \v e h\/(Jh- The degeneracy might be resolved by us- 
ing two beam lines with different 5*s. It is also possible 
to check the consistency of the observed DM wind with 
the astrophysical observation. 

In the above, we assumed that the beam energy is 
tuned to -Ebeam- The measurement of the selectron and 
neutralino masses with precision of the order of or be- 
yond an O(10) MeV level in the future collider experi- 
ments might be challenging, but important for this ex- 
periment. This can be seen in Fig. [21 where the modu- 
lation of the signal cross section is presented in cases of 
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FIG. 1: Modulation of cross section for e~x° — > e 
during a sidereal day for 8* = 0° and 45°. Here, we take 
m^o = lOOGeV, Am = lOGeV and Ebe&m = Bbcam- Other 
astrophysical parameters are given in text. 
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FIG. 2: Modulation of cross section for e x° ~~* e X° during 
a sidereal day for the beam energy -Ebeam = Bbcam — 20, 
-10, +0, +10, and +20MeV. Here we take m ? = lOOGeV, 
Am = lOGeV and 5* = 45°. 



£bcam = £bcam - 20, -10, +0, +10, and +20MeV. The 
cross section reduces significantly once \AE\ 3> 20 MeV 
for the parameters. If the error of the mass difference 
is more than 20 MeV, one has to scan the beam energy 
to find the signals. The measurement of the beam en- 
ergy dependence of the event rate might be useful to de- 
termine the -Ebcam, since the phases are reverse in the 
positive and negative energy deviation. Precise determi- 
nation of -Bbcam might allow us to interpret the event 
rate and determine the DM parameters. Especially, the 
measurement in the different beam energies may be used 
to resolve the parameter degeneracy between pdm and 
cr/j in the observed event rate. 

The SSIS model has two parameters, pdm and ah, in 
addition to the velocity and direction of the DM wind. 
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FIG. 3: Ratio of cross sections for e~ \ — > e~x in two cases 
that the beam axis is parallel and perpendicular to the DM 
wind, as a function of Am/m. Here, -Ebeam = Ebeam- 

We have seen that these might be constrained by chang- 
ing the beam axis and/or beam energy in the electron 
accelerator experiment, though a sizable event number 
is required. We stress that it might be possible when 
the selectron decay width is suppressed by the selectron 
and neutralino mass degeneracy or the coupling (or the 
neutralino mixing) and it is comparable to the typical 
deviation of y/s from the selectron mass in the DM ve- 
locity distribution. In Fig. [21 we show a ratio of the cross 
sections in two cases that the beam axis is parallel and 
perpendicular to the DM wind as a function of Am/m. 
Here, we use Eq. (j^J for the selectron decay width. Larger 
Am/m makes two cross sections closer since the selectron 
decay width is increased. In this case, the modulation of 
the cross section becomes featureless. 

The large-scale features of the flat rotation curves 
around galaxies are reproduced in the SSIS model. How- 
ever, numerous dynamical arguments suggest that actual 



halo model may not be well described by such a distri- 
bution. In addition to the axisymmetric |l7j ] and triaxial 
halo models 0], non-Maxwellian distributions, such as 
in the Sikivie caustic model 0, are proposed. The DM 
streams expected from the Sagittarius dwarf tidal stream 
might also affect the local DM velocity distribution [lfij . 
It is important to measure the local DM velocity distri- 
bution, including the directional dependence, so that the 
models are discriminated. Our proposal might be appli- 
cable to it. 

In this paper we discuss a possibility of neutralino DM 
direct detection through the scattering of local DM neu- 
tralinos with high intensity electron beam. The merit of 
the approach is that the DM can be identified as "the 
lightest neutralino" in direct and less ambiguous ways 
than the conventional DM detection experiments. Fur- 
thermore, the local DM density and velocity distribu- 
tion might be constrained. However, to study the neu- 
tralino DM, the current of the electron beam must be 
O(10) times higher than those currently planned, and 
the neutralino and selectron masses must be known very 
precisely at the level of precision expected in the future 
linear collider. In addition, the mass difference between 
them should be small to have a significant cross section. 
The experiment also requires the special setup to allow 
the long detector system along the beam line. 
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